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We show the existence of a pseudogap phase in the Nambu Jona-Lasinio model 
of quark interactions. In the pseudogap phase chiral symmetry is restored but qq 
pseudoscalar mesons still exist and they are massive. Such a behavior is interme- 
diate between a BCS superconductor and a Bose Einstein Condensate. We suggest 
the relevance of this phenomenon for an understanding of recent lattice QCD and 
experimental data. 
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I. INTRODUCTION 

Our current understanding of nuclear phenomena indicates that, at finite temperature, 
hadronic matter undergoes a phase transition to deconfined gluons and quarks. Quantum- 
Chromo-Dynamics (QCD) lattice calculations [1] strongly support the conclusion that, at 
some critical temperature, a transition to the quark-gluon plasma (QGP) phase occurs; this 
temperature is numerically equal to the critical temperature for the restoration of the chiral 
symmetry. 

The earliest suggested signature of QGP was the strong suppression of the light and 
heavy qq bound states for temperatures larger than critical temperature jy]. However, more 
recent lattice results 0, using the Maximal Entropy Method, have found that mesonic 
bound states, light or heavy, actually persist for temperature at least a factor of two larger 
than the critical temperature and the quasiparticle mass turns out to be still quite large |4| . 



Furthermore, the RHIC data on the radial and elliptic flows can be explained by partonic 
cascades j^| and viscosity corrections [?J only if the partonic cross section is about 50 times 
larger than the perturbative QCD calculations, which indicates a strong coupling regime. 

Some models have been proposed to describe this strongly interacting phase [a, Kg • The 
general feature resulting from both theory and experiment is the appearance of two different 
temperatures: T c and T* (T c < T*). The lower temperature T c should be associated with 
deconfinement and chiral symmetry restoration, the upper T* is related to the decoupling 
of the bound states from the spectrum [3j]. 

The presence of two temperature scales is an interesting phenomenon that has an ana- 
logue in high temperature superconductors. The nature of the phase transition for these 
systems is still matter of debate, see, for a recent review 10]. An established experimental 
fact seems however to be the existence in these superconductors of a pseudogap [Uj , which 
is a depletion of the single particle density of states around the Fermi level. Another char- 
acteristic feature of the high T c superconductors, is their coherence length £o which is much 
smaller than in ordinary superconductors. Both features might be the phenomenological 
manifestations of a crossover from the Bardeen-Cooper-Schrieffer (BCS) behavior of the or- 
dinary superconductors to a Bose Einstein Condensate (BEC) behavior [12]. In this scheme 
at a certain critical temperature a pseudogap phase is reached where the quasiparticles are 
still gapped although phase coherence and long range order are lost. Superconducting fluc- 
tuations might be responsible for this phenomenon, though other interpretations are also 
possible |1Q] . The relation between the pseudogap transition and field fluctuations has been 
qualitatively elucidated in the Nambu Jona-Lasinio (NJL) model jl^l . 14]; moreover it has 
been discussed in Il5j in the framework of a semi-phenomenological constituent quark model, 

n n 

in [16|] in the case of color superconductivity and in [17|] for low density nuclear matter. 

In this paper we analyze in more detail the dynamics of this NJL pseudogap phase, 
described by fermionic degrees of freedom and/or by the equivalent bosonic system, including 
also finite density effects. The picture that emerges from an explicit numerical calculation 
can be summarized as follows: 

i) Chiral symmetry is broken in the chiral limit, at T < T c and small density. The scalar 
meson, a, is massive and the pseudoscalar meson, tt, is the massless Nambu Goldstone Boson 
(NGB). 

ii) The temperature (and/or the density) restores the chiral symmetry and the order 
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parameter < ipip > goes to zero at T c . The signal of the restoration of the chiral symmetry, 
due to the background fluctuations, is given by the equal mass of the scalar and pseudoscalar 
mesons for T > T c . Therefore the NGB acquires a mass, while remaining in the spectrum. 
In this phase the constituent fermion mass remains finite (pseudogap regime) and there is a 
strong fermion-meson coupling. 

iii) Above another critical temperature, T* > T c , the constituent fermion mass goes to 
zero and the mesons decouple from the physical spectrum. 

The paper is organized as follows. In section |H] we present our formalism. In section ITTT1 
we discuss the bosonization of the NJL model and in section IIVI we compute the stiffness. 
In section El we compute the pseudogap and in section ED we draw our conclusions. 



II. CONSTITUENT QUARK MASS AND PION DECAY CONSTANT 

We use 

C = C + C 4 (1) 
with massless quarks at finite chemical potential in two flavors (Nf = 2). Here 

C = ^{id vl v + Wo)iJ> (2) 

while £4 gives the Nambu Jona-Lasinio (NJL) coupling (jl^ and for reviews 01): 

" C ''" UrW)* + (^75T^) 2 1 • (3) 



2N C 

In the sequel we neglect the so called exchange terms [3]. In the Mean Field Approximation 
(MFA) the constituent quark mass is obtained by the gap equation, which at T = fi = is 
as follows 

„ , T „ T G f d 3 p m g . t . 

m a = 4N f N c — / — V- ; q = , (4) 

f c 2AUa(2tt)3 



1 

where A is a 3D cutoff. One can note the relation between the constituent quark mass and 
the chiral condensate 

m, = -^(0|«u|0>. (5) 

( G \ 2 

Going beyond the MFA means considering diagrams of order I — — ) , see e.g. 
Clearly the smaller N c , the larger the role of the fluctuations. As a consequence, an effective 
momentum dependent quark mass is introduced M e ff(p). The momentum dependence varies 
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m,M eff (p) 



with the model, but for the model we use here, which corresponds to model A of 
is almost constant, with variations of a few percent in the whole p range. Therefore we will 
neglect these effects altogether. Moreover it can be observed that the reliability of the 
mean field gap equation to evaluate the constituent quark mass also in a region where the 
fluctuations are relevant (as for example for small N c ), has been explained by Witten (2^ 
in exactly solvable models in 1+1 dimension and numerically checked by lattice calculations 
in 2+1 dimensional models j^J. 

The gap equation at finite T and /1 can be obtained by considering a sum over Matsubara 
frequencies, by the substitution of the energy integration variable p° with iu n — /i, where fi 
is the quark chemical potential and u n = nT(2n + 1). 

Eq. (0J) relates the constituent quark mass m q at T = fi = to the NJL coupling Go 
and the cutoff A. Its generalization for finite T and fi provides the T and //-dependent 
constituent quark mass in the mean field approximation: m q (T,fi). By getting rid of the 
NJL coupling constant, we get m q (T, fi) from the constituent mass at T = fi = as follows: 







d 3 p 

A (2+)3 



sinh y 



^jp 2 + m\ e (cosh y + cosh x) 



(6) 



where 



fi e y / p 2 + m2(T,/i) 
z = 7^ , V = 7^ , e = tf, • ( 7 ) 



T T 
The pion decay constant at T = \i = is 

fl = -Um\N c f l - Y — 2 . (8) 

y ja (27r) 4 (p 2 — 171^ + ie) z 

One can use this equation to get A by fixing the constituent quark mass. Using m q = 300 
MeV as an input, from f n = 93.3 MeV one gets A = 675 MeV. 

We define the T* = T*(/x) temperature by m q (T*,fi) = 0; numerically one finds, at 
fi = 0, T* = 185 MeV. For generic values of \x the result can be obtained by Eq. ((HJ). The 
generalization of Eq. (jHJ) to finite T and \x is straightforward. 



III. BOSONIZATION OF THE NJL 



Let us first work at T = fi = 0. As is well known the NJL model can be made equivalent 
to the linear a model 24| (see also Q]). One introduces fields a and 7r by writing the 
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generating functional as follows 



Z = I [#][#] exp \ I dx(C + £ 4 ) \ = I [#][#] [dcr)[d7r l 

9 2 N C 



x exp || dx U(id v -f + /ijo ~ 9o{v + HnT ■ - ^ [a 2 + tt 2 ] J 1 (9) 



where 



9o = ^ (10) 

is the meson-quark coupling constant. Eq. (fTU|) is the analogous of the Goldberger-Treiman 
relation. By bosonization and derivative expansion the NJL model becomes equivalent to 
the a model with lagrangian 

C = £ ( (daf + (d*) 2 + £ (a 2 + tt 2 - M . (11) 



2 \/ ' v ' 4 
The parameter /3 is given by j^J, 

/3 = % 2 iV c iV / (/ -2fi ) , (12) 

with 

d A p E 1 



(2tt)* (p| + 



1 /• rf 4 p g p| 
"° " 4 7 (2tt)Mp| + m2)3 • ^ 

Classically one has < 7r a >= 0, < a >= f n ; as a consequence the field a' = a — f n acquires 
a mass nf n / a/2 while the pions remain massless. 

If the a' mass is very large, the model is equivalent to the non linear a model with action 

S= l ljd A x ((da) 2 + («9tt) 2 ) (14) 

and fields satisfying the constraint 

- 2 + - 2 = /. 2 - (15) 



IV. STIFFNESS 



The generating functional of the nonlinear a model for vanishing external sources can be 
written as follows 

Z= [[da][diT}[d\}exp{iS} (16) 



6 



S = ^Jd D x ((da) 2 + (d*) 2 — A (a 2 + 7r 2 — f 2 )) (17) 

where the A functional integration implements the non linear condition on the fields (j!5|) . 
After integration over the iz fields we get 



with 



[da] [dX] exp(iS e ff[a, A]) 



S eff [a, A] = *-Tr \n[d 2 + X(x)) + t J ' d D x (f 2 X(x) - a[d 2 + X(x)]a 



(18) 



(19) 



If the number of flavors is large we can use the saddle point approximation and search 
for solutions a, X independent of x. Since we are looking for phase transitions at finite T 
and /i, we consider Matsubara frequencies and non vanishing baryonic chemical potential. 
One gets in this way 2^|: 

d 3 k T 



= Act 



(2vr) 3 {u n - ifx) 2 + k 2 + A 



(20) 



with iu n = 2nnT; therefore one has two solutions 1 25]: 



i) a = and A implicitly given by 



ii) A = and 



oo „ 

Pf* = (N 2 -l) £ / 



N 2 - 1 00 

<? = K-- L ir- £ 



d 3 k 



T 



(2vr) 3 (£o n - z/i) 2 + k 2 + X 



d 3 k 



T 



(21) 



(22) 



(3 J^J {2<KY{u n -in) 2 + k 2 

The case i) corresponds to chiral symmetry restoration, the case ii) to spontaneous symmetry 
breaking and massless pions. Note that in the chiral symmetric phase, a = and also 
(0\qq\0) = 0. 

Let us now consider the case of finite T and fi. The critical line of chiral symmetry 
restoration is defined by A = o = 0, or 

d 3 k T 



Nf-1 ~ 
f 2 (T, A*) 8tt 2 Jo 



(27T) 3 (cD n -2/i) 2 + fc 2 

' fi + k\ 



dk k 



coth 



2T 



coth 



' fi — k" 
2T 



(23) 
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At finite temperature and density we have for the stiffness, from Eqns. (jl2j) . ()13|) 

d 3 p 



P = 2g 2 N f N c T J 



2 -L- (w„ - ^) 2 ] 2 [e 2 + (w n - z/x) 2 p 



(24) 



/a (2tt) j _ 

The critical temperature T c as a function of [l is obtained by equating (}2*3*j) and (|24|) jl^ . 
jl^ . The value of T c depends on the value of the pionic cutoff A^; typical values should be 
of a few units of f n since f n fixes the scale of the derivative expansion. Numerically we find, 
at n = 0, for K = 200 MeV, T c = 146 MeV; for = 300 MeV, T c = 127 MeV. We stress 
the assumption of two different cutoffs, A and A^; for a discussion on this point see I VII 

If T c (/x) defines the critical line, we can compute the mass of the composite boson formed 
by the qq pair above T c (/z), i.e. in the Wigner phase. The existence of such pairs in the NJL 



model is a feature that was observed long ago [26j and interpreted as precursor of chiral 
phase transition at finite temperature and chemical potential. We differ from these authors 
for two important aspects. First, we identify the critical temperature as T c and not as T*. 
Second, in 2(| the fluctuations are identified by looking at peaks in the u distribution of 
the strength function 

A(k, uj) oc ImG R (k, to) (25) 

at k = (G R is the retarded Green function describing the fluctuation of the order parameter 
in the Wigner phase). Here we use the stiffness to take into account the fluctuations and 
compute the mass of the ir mode m 2 = A by identifying /3(A) in (}2*T]) with (3 given by ()24|) . 

For n = the results are reported in fig. ^for the value of the parameter A^ = 200 MeV 
and for N c = 3 (solid lines) and A^ c = 10 (dashed lines). Before discussing the dependence 
on A^ c , let us note that in fig. ^the function er(T) differs from zero in the interval (0,T C ) 
and vanishes in (T c , T*) for both values of N c . On the contrary the function m T (T) vanishes 
in (0, T c ) and rapidly increases in (T C ,T*) with a divergent behavior around T*. 

An interesting aspect of these results is that they allow to define a cross-over temperature 
Tpair > T c ; in our computation T pair is the temperature at which the composite bosons 
disappear from the physical spectrum. As a result of our approximation we find T pa i r ~ T*. 
Lattice QCD investigations show that the deconfinement and chiral phase transitions take 
place at the same temperature, see e.g. and H- 

Though our results indicate the 
existence of two different temperatures, they do not contradict lattice studies. As a matter 
of fact, the NJL model does not incorporate confinement. Therefore we can only state that 
our lower critical temperature T c is the chiral restoration critical temperature, while leaving 
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FIG. 1: - The parameter a and mass of the pionic mode as functions of the temperature T 
at li = (in MeV). T c is the critical temperature. For T > T c a vanishes, while vanishes for 
T < T c . Values of the parameters are = 200 Me V, A = 675 Me V. Solid lines refer to the case 
N c = 3, dashed lines to N c = 10 (weaker coupling). As a result of the computation, T c ~ 146 MeV 
for N c = 3 and T c ~ 172 MeV for N c = 10. T* ~ T pa i r is the dissociation temperature of the qq 
pairs. 

undecided if it coincides with the deconfinement temperature, defined for example as the 
temperature where the string tension vanishes. As for the upper characteristic temperature 
T pair ~ T*, it is unrelated with the deconfinement temperature and is defined either as 
the temperature where the strongly correlated Cooper pairs, behaving as massive mesonic 
states in the interval (T c , T*), decouple from the spectrum, or as the temperature where the 
pseudogap phase disappears. 

We conclude that the behavior of this model is more akin to a Bose-Einstein condensate 
than to a BCS system. Indeed while in the latter case the temperature characterizing pair 
formation coincides with the condensation temperature, the former case is generally charac- 
terized by two critical temperatures, one (upper) temperature when the pairs are preformed, 
and another (lower) where the pairs condense (for discussions in condensed matter see e.g. 

Q n 
and for a review [10]). This behavior is confirmed by the large N c limit of the present 

model. In this limit the four fermion coupling becomes weak, the mean field approximation 

is exact and the two temperatures merge. This behavior can be easily proved by taking 



9 



the large N c limit of the equation for the stiffness that determines T pair . For N c — > oo one 
gets T c ~ T pair , and the model defines a BCS superconductor. By way of example we have 
reported in Fig. ^the result of the calculations for N c = 10 (dashed lines) showing that for 
weak coupling T c — > T*. 



V. PSEUDOGAP 



The pseudogap can be seen by the density distribution N(cu): 

dk 



N(u) = 4j—Tr cJ p\k,u) (26) 



(c, / =color, flavor), where 
(d =Dirac) and 

G (k, lo) is the analytic continuation of the imaginary time Green function Q, in other words 



p°(k,u) = ^Tr dl A(k,u) (27) 



A(k, u) = —ImG R (k, u) . (28) 

7T 



G R (k,u;)=[G- 1 (k,uj + ie)-E R (k,uj)] 1 (29) 

where G~ 1 (fc, u) = (u + Ai)7° — k ■ 7 — m q is the analytic continuation of the imaginary time 
Green function for free quarks Q~ l {k, to) = (iu n + fi)-y° — k ■ 7 — m q . Notice that including 
only the mass term, i.e. at the lowest order in the interactions, we have 

iV Nt 1 

N(u) = N (u) = + ^{u + nY-m* . (30) 

We include the quark mass term in G(k,u) in the spirit of an effective lagrangian approach, 
where the four-fermion coupling modifies the free quark dispersion law and the interactions 
are provided by the fluctuating pion field. Actually the absence of interactions would imply 
also the vanishing of m q and one would get, in this case, 



N(u) = N free {uj) = -^{uj + fx) 2 . (31) 

To compute the effect of the interactions at the next order we consider the effective lagrangian 
displayed in eq. (jOJ): 

C eff = $[id u Y + fi l0 -g (<T + ilr>r ■ tt)]V - 7^ k 2 + 7T 2 1 (32) 
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together with the condition (|15|). At the lowest order in the pion field we have the following 
interaction lagrangian 

tint = $ ^^J" 71 " 2 ~ ^075^ • T^j 4> (33) 

together with the mass term —g^j^^ = —m q ipip. At the lowest order in g$ we find 

E R (k, uj) = i ((«! + ia 2 ) lo + (Ai + iA 2 )k • 7 + Ci + Kz) (34) 
and (ai, Aj , Ci real): 

a 1 + ia 2 = -(JV?-1)T g / [«(q* " " " ^ 

n=-oo 

Ai + ^a 2 = -(jv?-i)r £ / [*]flg(q 2 , n 

n=-oo JA7r ^ 

Ci+< 2 = (iVj-l)T £ / [dq}gl{^ [ 1 m 



' Aw d 

1 - q • k/k 2 



^ 

2m g D 



(35) 



The measure [dg] is 

d 3 g 1 
= (27T)3(27™T) 2 + q 2 + m2 (36) 

and 

D= [u n + i(u + n)} 2 + (k-q) 2 + m 2 q (37) 
with o; n = lisTin + 1/2). We have included a form factor in the definition of g^: 



^'> = * U + (Wq' J ' <38) 
It corrects the large g 2 behavior of the loop diagrams and plays a role analogous to the 
regulator used in |25j. We take mp = 100 MeV, but the exact numerical value is of no 
interest here as we are only interested in the qualitative behavior of the density of states. 
We note that the effect of the form factor is to reduce the perturbative contribution to 
the density distribution N(u). As a matter of fact this perturbative term decreases with 
decreasing mp, while for nip — > 00 the effect of the form factor vanishes. Our numerical 
choice, mp = 100 MeV, only indicates an order of magnitude. Much larger values are 
excluded since they would extend the validity of the approximation beyond the kinematical 
range dictated by the neglect of higher order terms in the derivative expansion of eq. (jllj) . 
We also note that a further effect of the form factor is to make the series in eq. ()35|) 
rapidly convergent, which allows an inversion between analytic continuation and Matsubara 
summation, see e.g. ref. and related work in refs. I30I \&\. 
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FIG. 2: - The ratio N(u)/N f ree (to) as a function of the energy to. N(uS) is the density of states for 

the interacting theory; Nf ree (to) the density for the free fields. Results are presented for [x = 10 MeV 

and for various temperatures. The solid line represents the absence of interactions and corresponds 

T — T c 

to T = T* . The other curves refer to three different values of 5 = — — — : 5 = 0.02, dotted; 
5 = 0.15, dashed and 5 = 0.20, dotted-dashed line. 



Similarly to the Nf ree case, the |k| integration can be performed using the residues the- 
orem. Our results are as follows; we define N(u) = N (u) + N pert (uj) and we get 

N pert {uj) = — ^- k a 2 {uj, ho) H — 1 . (39) 

7T Z V ZKq Z J 

Here /' denotes partial derivation in the k variable, 



k = y/(u + tA)*-m? q (40) 

and 

/(w, k) = -2(w + fi)[(u + ti)a 2 (w, k) + k 2 X 2 (u, k) + m q ( 2 (uj, k)} . (41) 

Our results are displayed in Fig. |21 where we plot y for three different values of the 

Nf ree (u) 

ratio 

T — T 

-L c 

i.e. 5 = 0.02 (dotted line) 5 = 0.15 (dashed line) and 5 = 0.20 (dotted-dashed line) and 
for one value of the chemical potential (fi = 10 MeV). For other values of \x the results 
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are similar. The figure shows the existence of a pseudogap in the (T c , T*) region. The 
pseudogap is maximal near T c and decreases for T c —>■ T*, which can be easily understood 
as the combined effect of the vanishing constituent mass m q and the divergent behavior of 
run for T — > T*. 



VI. CONCLUSIONS AND OUTLOOK 



The previous analysis shows that in the T—fi phase diagram of the NJL model, at a fixed 
value of \x (smaller than the critical value when the effects of color superconductivity become 
significant), there are two typical temperatures, associated with chiral symmetry restoration 
and with the dissociation of the qq pairs. This qualitative behavior is suggested in ^| by 
Babaev, who changes the conclusions of 13( employing a similar method. The authors of 
ref. 1 131 use A w & A and conclude that the NJL model has no chiral broken phase, unless the 
number of colors N c is larger than its physical value N c = 3. We agree with |l4j where two 
different ultraviolet cutoffs are used, as they are related to different degrees of freedom. Our 
explicit calculations corroborate therefore the conclusions of [ljj . In particular we explicitly 
prove the appearance of a pseudogap, related to the existence of a strong coupling. The 
existence of two different cutoffs, one fermionic (A) and another one bosonic (A w ) is due to 
the fact that, as shown in refs. Q and Q, the regularization of quark loops by means of A 
does not restrict the four momenta of the mesons. In all non-renormalizable field theories, 
such as the NJL model, new parameters are introduced as more loops are included in the 
calculations. Since the inclusion of internal boson lines corresponds to add more loops, the 
presence of meson propagators implies the existence of the new cutoff A„. 32|. In 33j the 
interval (0, 1.5) for A n /A is investigated, but in (see alsoj^l) it is shown that > A 
leads to instabilities. Our choice A„- < A is consistent with 34]: in particular we find that 
the value A^ =200 MeV produces a numerical value for the chiral phase transition in rough 
agreement with lattice results. We finally note that the existence of two different cutoffs can 
be also interpreted as the consequence of the non-universality of a critical stiffness in the 
non linear sigma model for 3+1 dimensions. This lack of universality implies 3] that the 
knowledge of (3 C in the NJL model does not allow, by alone, to fix the position of the phase 
transition in the effective field theory. Therefore the appearance of a new parameter appears 
quite natural from this viewpoint. Summarizing: Our derivation is based, differently from 
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Ref. Q], on two cutoffs and, as such, it mi ght be controversial. However, on the basis of 

14 1 summarized above, we believe that using 



the arguments of Refs. 



32 



33, 



34 l35| and 



two different cutoffs is a proper procedure. 

As already mentioned in the introduction, some other models have been proposed to 
understand in QCD the persistence of the qq bound states above the critical temperature and 
the hydrodynamical properties of the medium observed at RHIC; let us briefly comment on 
them. In j8j a strong QCD Coulomb-like interaction and the formation of a multiple bound 
states of quasiparticles has been suggested as a possible explanation of the data. However 
it is not clear if this stron g C oulomb-like interaction is obtained in lattice simulations above 
the critical temperature jflfij ]. More recently in j^] a strong non-confining qq potential is 
predicted in the framework of the field correlator method, which might account in part for 
the data. 

The NJL model is obviously different from QCD; in particular it does not contain one of 
the key dynamical QCD effects, i.e. confinement. Nevertheless it contains the most relevant 
dynamics of the effective four fermion interactions and is largely used as an effective model of 
low energy QCD. Therefore the similarity between its behavior and some non-perturbative 
features of QCD (lattice results and the phenomenological interpretation of the RHIC data) 
should not cause much surprise. On the basis of the present work, it would be interesting to 
investigate if a pseudogap also shows up in lattice QCD at intermediate temperatures and 
small chemical potential. 
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